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When  a  number  of  TEMs  (thermoelectric  modules)  are  connected  in  a  series— parallel  matrix  and  under 
mismatched  temperature  gradients,  the  overall  maximum  output  power  of  the  thermoelectric  generator 
(TEG)  may  be  lowered  by  certain  TEMs  with  relatively  smaller  temperature  difference.  It  is  possible  to 
avoid  such  a  performance  decrease  by  the  disconnection  of  these  low  temperature  TEMs,  provided  that 
the  critical  power  point  can  be  accurately  determined.  In  this  paper,  firstly  a  rigorous  and  universal 
formulation  is  fully  detailed  to  mathematically  determine  the  conceptions  and  conditions  of  the  critical 
power  point  in  the  series  and  parallel  TEM  arrays.  Secondly,  experiments  of  a  series— parallel  hybrid 
interconnected  TEG  are  presented  to  clearly  quantify  the  theoretical  analyses.  Finally,  the  hierarchical 
simulation,  based  on  the  SPICE  (simulation  program  with  integrated  circuit  emphasis)  platform,  is 
applied  to  estimate  the  critical  power  point.  By  numerically  modeling  the  nonlinear  physical  processes  of 
the  TEG,  the  simulation  can  be  used  as  an  enabling  technique  in  any  model-based  controller  to 
dynamically  minimize  the  mismatch  power  loss  within  the  TEM  matrix  of  any  configuration.  In  exper¬ 
imental  and  numerical  results,  a  number  of  critical  power  points  are  disclosed  for  a  2  x  4  parallel— serial 
hybrid  TEM  matrix,  where  the  hot  temperature  mostly  ranges  from  120  °C  to  60  °C. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

TEGs  (Semiconductor  thermoelectric  generators)  have  the 
merits  of  high  durability,  great  environmental  harmony,  and  pure 
DC  (direct  current)  power  source.  As  they  can  recover  the  huge 
amount  of  low  grade  waste  heat  into  electricity  in  a  simple  manner, 
TEG  systems  have  been  of  great  interest  to  energy  applications  in 
recent  years.  Representative  examples  of  TEGs  include  fuel  fired 
combustors  and  thermal  fluid  tubes  in  industrial  power  plants  and 
private  sectors  [1-7],  as  well  as  the  exhaust  pipes  in  fuel  cells 
[8-11],  aircrafts  [12]  and  vehicles  13-19],  where  thermoelectric 
modules  (TEMs)  can  be  mounted  on  heat  sources  to  operate  as  an 
array  of  thermal  batteries  that  supplies  power  to  the  load. 

In  composing  a  battery  pack,  ordinary  voltage  sources  such  as 
alkaline  batteries  usually  prefer  the  series  connection.  This  is 
because,  if  the  batteries  are  mismatched,  then  those  with  higher 
voltage  may  charge  others  when  the  batteries  are  connected  in 
parallel.  In  this  case,  the  batteries  with  lower  voltage  cannot  work 
normally  to  output  power  and  may  be  overheated  by  the  great 
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current  as  the  resistance  of  most  batteries  is  very  small.  The  ther¬ 
mal  runaway  may  even  be  aroused  for  lithium-ion  batteries  by 
short  circuit,  that  both  heat  and  temperature  are  kept  increased  to 
ultimately  lead  battery  packs  to  explode.  In  general,  a  lot  of  voltage 
sources,  e.g.,  battery  charger,  backup  power  source,  and  grid- 
connected  generators,  can  only  be  paralleled  after  special  protec¬ 
tion  design  to  prevent  such  mismatch  damages.  Different  from 
these  voltage  sources,  mismatched  TEMs  do  not  have  the  charging 
and  security  risks  in  the  case  of  parallel  connection.  Consequently, 
series  and  parallel  connection  are  of  equal  importance  in 
composing  TEM  arrays  to  satisfy  the  technical  tradeoff  between  the 
reliability  and  the  power  generation  specifications  imposed  by  the 
practical  load  voltage  and  current. 

Yet  a  well  recognized  rule  in  designing  such  a  series/parallel 
hybrid  battery  pack  as  TEM  arrays  is  that  every  cell  should  have 
close  electric  parameters  in  order  to  make  the  power  loss  caused  by 
the  mismatch  minimal  [20].  The  electric  parameters  refer  to  the 
open-circuit  voltage  and  the  internal  resistance  for  Thevenin's 
voltage  sources.  As  the  TEMs  connected  together  are  usually  rep¬ 
resented  by  Thevenin's  voltage  source  model,  they  should  ideally 
be  of  similar  voltage  and  resistance  characteristics  to  each  other  so 
that  no  degraded  battery  limits  the  string  current  or  the  total 
voltage.  However,  both  theoretical  [1-4,6,7,9,11,12,16,18]  and 
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experimental  [5,8,10,13-15,17,19]  investigations  clearly  show  that 
the  variance  of  the  temperature  difference  applying  to  different 
TEMs  in  most  thermal  systems  can  be  as  much  as  a  few  hundred 
degrees.  Since  the  voltage  as  well  as  the  internal  resistance  of  a  TEM 
are  nonlinearly  determined  by  its  hot  and  cold  temperatures,  the 
TEM  mismatch  power  loss  is  inherently  associated  with  most  TEM 
arrays. 

As  can  be  proven  in  terms  of  the  fundamental  circuit  theory  and 
shown  by  means  of  numerical  calculation  [3,7,11]  and  experiment 
test  [21  ],  the  maximum  power  capability  of  a  TEM  array,  measured 
by  the  ratio  between  its  equivalent  open-circuit  voltage  and 
equivalent  internal  resistance,  may  be  lowered  by  the  inclusion  of 
some  TEMs  with  smaller  temperature  difference  due  to  the 
mismatch  power  loss.  In  this  case,  however,  the  possibility  to 
optimize  the  array  output  power  by  the  removal  of  the  TEMs 
operating  with  lower  heat  source  temperature,  implemented  by 
replacing  the  traditionally  fixed  hardwire  interconnections  be¬ 
tween  the  TEMs  with  controllable  switches,  has  only  received  some 
attention  recently  [22].  Once  the  array  is  spitted  into  two  sub¬ 
arrays,  and  one  can  be  assumed  as  the  original  array  whilst  the 
other  is  composed  of  the  removed  TEMs,  it  is  also  possible  to 
configure  each  of  the  sub-arrays  with  an  individual  power  elec¬ 
tronics  stage  and  MPPT  (maximum  power  point  tracking) 
controller  [23].  In  this  way,  not  only  the  mismatch  power  loss 
caused  by  the  temperature  nonuniformity  is  recovered,  but  more 
power  may  also  be  produced  than  the  original  TEM  array  with 
uniform  temperature  distribution  and  zero  power  loss. 

From  the  practical  viewpoint,  it  is  thus  important  to  both 
qualitatively  and  quantitatively  evaluate  how  small  the  tempera¬ 
ture  difference  of  certain  TEMs  is  enough  that  their  disconnection 
from  the  array  would  be  beneficial.  The  issue  has  not  been  thor¬ 
oughly  covered  by  the  previous  works  focusing  on  the  controller 
design  and  implementation  22,23],  but  can  be  identified  in  an 
alternative  way  to  find  out  the  extremum  of  the  maximum  output 
power  for  a  TEM  array  with  various  temperature  distribution  pat¬ 
terns.  Such  a  power  maximum  determines  the  best  occasion  to 
optimally  remove  TEMs/split  the  array,  and  is  corresponding  to  the 
so-defined  critical  power  point.  The  main  objective  of  this  paper  is 
to  carry  out  the  qualitative  and  quantitative  evaluation  mentioned 
above,  and  the  diagnoses  of  the  critical  power  point  will  be  realized 
by  means  of  the  analytical,  experimental,  and  numerical  ap¬ 
proaches,  respectively,  aiming  for  the  feasible  design  assistance  and 
guide  in  building  practical  TEG  systems.  The  main  contribution  of 
the  analytical  study  is  to  rigorously  analyze  the  mismatch  power 
loss  in  TEGs  and  formulate  the  critical  power  point  to  optimally 
split  TEM  arrays  in  terms  of  both  series  connection  and  parallel 
connection  (Section  2).  The  main  contribution  of  the  experiments, 
taken  from  a  practical  TEG  consisting  of  eight  independently 
controlled  TEMs,  is  to  confirm  the  theories  of  critical  power  point 
and  provide  concrete  figures  on  the  magnitude  of  thermal  imbal¬ 
ance  to  arise  such  critical  power  points  (Section  3).  The  main 
contribution  of  the  thermo-electro  coupled  simulation,  using  the 
hierarchical  TEG  modeling  in  the  circuit  simulator  SPICE  as  the  ideal 
solution,  is  to  accurately  predict  TEG  behaviors  at  the  system  level, 
where  multiple  TEMs  can  be  both  electrically  connected  and  ther¬ 
mally  combined  arbitrarily  (Section  4). 

2.  Mismatch  power  loss  and  critical  power  point 

A  common  configuration  to  create  the  temperature  difference 
for  a  TEG  is  the  fluid-to-fluid  heat  transfer,  where  a  number  of  TEMs 
are  placed  in  the  heat  exchanger  interface  separating  a  hot  fluid 
stream  and  a  cold  fluid  stream,  and  generate  power  as  the  heat 
flows  through  them  [1-4,7,10-14,19].  Each  TEM  is  electrically 
described  as  a  voltage  source  serially  connected  with  a  resistor  [24], 


and  when  all  TEMs  are  connected  in  series,  the  total  voltage  and 
resistance  of  the  TEG  in  the  quasi-one  dimensional  (ID)  configu¬ 
ration  [25]  can  be  described  as. 
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As  shown  in  Eqs.  (1)  and  (2),  the  local  power  generation  of  any 
TEM  in  the  TEG  system  can  be  developed  by  the  nonlinear  integral 
functions  within  the  local  temperature  difference,  i.e.,  from  Tc(i)/ 
0 (i)  to  Th(i)/L(i),  where  heat  source  temperatures  Th  and  Tc  can  be 
assumed  same  for  all  thermocouples  of  each  individual  TEM.  The 
power  capability  of  the  TEG  with  all  TEMs  connected  together  is 
described  as. 


P  V2*1 

(R  +  Ri)2’ 

which  is  maximized  when  R\  =  R ,  i.e., 


(3) 


V2 

Pmax  =  4^-  (4) 

Similarly,  when  all  TEMs  are  connected  in  parallel,  V  and  R  of  the 
TEG  can  be  described  as. 


According  to  Kirchhoff  circuit  laws,  the  numerator  of  Eq.  (5) 
represents  the  sum  of  the  ratios  between  the  voltage  and  the 
resistance  of  all  TEMs,  i.e.,  XI  and  the  denominator  of  Eqs.  (5) 
and  (6)  represents  the  sum  bf  the  reciprocals  of  the  resistance  of  all 
TEMs,  i.e.,  X 1  /£/•  Although  all  TEMs  are  connected  in  parallel,  the 
N  pairs  of  thermoelements  within  each  TEM  are  usually  connected 
in  series,  as  shown  in  (5)  and  (6). 

Other  instances  of  TEG  heat  source  include  various  waste-heat 
recovery  applications  in  the  outer  walls  of  exhaust  pipes  or  com¬ 
bustors  [5,6,8,15-18].  When  all  TEMs  are  connected  in  series,  the 
total  voltage  and  resistance  of  the  TEG  in  the  two-dimensional  (2D) 
configuration  can  be  similarly  described  as. 
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Eqs.  (3)  and  (4)  can  still  be  applied  to  obtain  the  power  output  in 
the  2D  case. 

The  circuit  of  the  series  connected  TEMs,  described  by  Eqs.  (1), 
(2),  (7)  and  (8),  can  also  be  generally  visualized  by  Fig.  1  (a),  in 
which  M  may  be  equal  to  |lf/w|  or  [Dw/S\  for  the  ID  or  the  2D 
case,  respectively.  The  voltage  source  Vj  composed  by  the  dotted 
line  represents  the  depressed  TEM  subjected  to  a  small  A 7}.  If  Vj  is 
idealized  that  its  resistance  Rj  =  0,  the  power  output  of  the 
remaining  TEMs  in  the  circuit  system  will  not  be  affected  irre¬ 
spective  of  how  small  Vj  is.  Nevertheless,  practical  TEMs  always 
exhibit  an  internal  resistance  that  in  this  way  Rj  will  lower  the 
power  capability  of  other  high  voltage  TEMs,  and  may  even  cause 
the  TEM  string  with  Vj  to  deliver  less  power  to  the  load  than  the 
remaining  TEMs  without  Vj. 

In  order  to  investigate  the  condition  for  such  a  critical  power 
point,  e.g.,  how  small  Vj  or  how  big  Rj  is,  it  is  necessary  to  rigorously 
formulate  the  power  loss  caused  by  the  mismatch  among  the  TEMs 
connected  together.  By  reducing  Eqs.  (1),  (2),  (7),  (8)  and 


(b) 

Fig.  1.  Power  source  pack  with  TEM  cells,  (a)  Series  connection  model,  (b)  Parallel 
connection  model. 


substituting  them  into  Eq.  (4),  the  comparison  between  the  TEG 
Pmax  and  the  sum  of  Pmax  of  all  TEMs  yields. 


According  to  Cauchy  inequality,  Eq.  (11)  is  equal  to  zero  only  if 
Vi/fti  =  V2IR2  =  ...  =  Vj/Rj  =  ...  =  VmIRm,  which  means  that  the  TEG 
Pmax  is  equal  to  the  sum  of  Pmax  of  all  TEMs  when  they  generate 
individually;  otherwise  Pmax  <  )Ci=i Pmax  v  and  the  difference  be¬ 
tween  them  is  the  mismatch  power  loss  occurring  inside  the  TEG. 
Considering  the  simplest  circuit  of  two  TEMs,  if  both  V\  ¥=  V2  and 
Pi  =£  P2,  in  principle  V1/P1  =  V2/P2  is  still  possible  that  the  equality 
holds  in  Eq.  (11).  However,  according  to  Eqs.  (1),  (2),  (7)  and  (8)  as 
well  as  the  temperature  dependencies  of  aPtYl  and  pp  n  of  Bismuth 
Telluride  26],  the  mostly  commercialized  thermoelectric  material, 
Pi  2  will  increase  when  the  AT12  becomes  lower  but  Vjt2  will 
decrease  for  the  same  change.  The  condition  V1/P1  =  V2/P2  to  allow 
such  a  lossless  mismatch  is  consequently  of  little  probability  to 
appear,  not  to  say  the  practical  series  connection  with  more  TEMs. 
Thus,  a  TEG  without  any  power  loss  usually  consists  of  matched 
TEMs  with  the  same  V  and  P  due  to  the  same  AT;  if  the  TEMs  in  a 
series  connection  are  mismatched  due  to  the  nonuniform  tem¬ 
perature  distribution,  there  is  usually  a  power  loss. 

If  the  TEG  in  Fig.  1(a)  cuts  off  the  low  voltage  TEM  Vj  and  other 
cells  in  the  TEM  string  operate  with  identical  AT,  then  Pmax  becomes 
equal  to  the  power  sum  of  all  the  remaining  individual  TEMs 
because  there  is  no  mismatch  power  loss.  The  condition  for  the 
critical  power  point  can  be  written  as. 


where  the  terms  within  the  parenthesis  amounts  to  the  new  Pmax 
without  Vj.  Eq.  (12)  can  be  more  clearly  represented  as 


f  ‘ £  (£"'< '  vf 

h4R<  4£",B,  ~4Ri 


(13) 


where  the  left  side  is  simply  the  mismatch  power  loss  introduced 
by  Eq.  (11).  Therefore,  when  TEM  Vj  is  shorted,  although  the  TEG 
terminal  voltage  would  be  lowered,  the  thermal  battery  pack  may 
perform  better  as  the  power  loss  caused  by  the  weakest  cell  in  the 
pack  is  greater  than  its  power  contribution. 

Similarly,  the  circuit  of  TEMs  connected  in  parallel  as  well  as  the 
depressed  TEM  Vj  due  to  the  non-uniform  temperature  distribution 
is  shown  in  Fig.  1(b).  The  identical  end  voltage  of  the  paralleled 
TEMs  and  their  built-in  Seebeck-voltage  imbalance  will  lead  to  a 
certain  amount  of  power  loss,  that  the  TEMs  with  different  po¬ 
tentials  waste  a  part  of  their  work  to  counteract  against  each  other. 
The  counteracting  phenomenon  is  particularly  obvious  when  the 
load  is  in  open  circuit  and  no  power  can  be  extracted,  where  there 
will  still  be  a  loop  current  internally  flowing  through  the  TEMs,  as 
shown  in  Fig.  1(b).  Because  Vj  is  less  than  the  voltage  of  other  TEMs, 
the  loop  current  flows  from  its  anode  to  its  cathode,  but  flows 
oppositely  from  the  cathodes  to  the  anodes  in  the  remaining  TEMs. 
In  the  same  way,  the  branches  of  each  paralleled  segment  will  have 
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different  currents  in  the  case  of  a  closed  load.  Such  phenomena 
within  the  power  source  caused  by  the  AT  mismatch  among  the 
parallel  TEMs  will  decrease  the  power  performance,  as  the  TEM  Vj 
partly  performs  as  a  load  for  the  power  source,  decreasing  the 
loading  capability  of  the  overall  TEG. 

The  power  loss  caused  by  the  mismatch  among  the  TEMs  con¬ 
nected  in  parallel  can  be  obtained  by  reducing  Eqs.  (5),  (6),  (9),  (10) 
and  substituting  them  into  Eq.  (4). 


M 


P max  —  P] max  i  — 

i=  1 


\  ~\M  V[ 

^'=lfy 
4  YM  - 
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vt 


--Y— 

^4 Rj 
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fv-M  _Vi _ 2_\ 

‘  (J_)2  "Mw  - 


(14) 


By  applying  the  Cauchy  inequality  theory  again,  Eq.  (14)  is  equal 
to  zero  only  if  V\  =  V2  =  ...  =Vj=  ...  =  VM,  where  the  TEG  Pmax  is 
equal  to  the  sum  of  the  individual  Pmax  of  all  TEMs  and  there  is  no 
mismatch  power  loss;  otherwise  Pmax  < Erli^max  v  and  the  dif¬ 
ference  between  them  is  the  mismatch  power  loss  occurring  inside 
the  TEG.  If  the  TEG  in  Fig.  1(b)  cuts  off  the  TEM  Vj  and  other  TEMs 
operate  with  identical  AT,  Pmax  becomes  equal  to  the  power  sum  of 
all  the  remaining  individual  TEMs.  The  condition  for  the  critical 
power  point  of  the  parallelly  connected  TEMs  can  be  written  as. 


>0, 


(15) 


where  the  terms  within  the  parenthesis  amounts  to  the  new  Pmax 
without  Vj.  In  terms  of  Vj  and  Rj,  the  critical  condition  can  also  be 
written  as 


M  V2  Eill£  W 

h  4Ri  -4  Rj  J 

where  the  left  side  is  the  mismatch  power  loss  when  the  branch  of 
TEM  Vj  is  not  open. 

The  previous  works  [3,7,11,21  ]  have  validated  the  critical  power 
point  of  the  TEMs  connected  in  series  by  means  of  numerical 
calculation  and  experiment  test.  In  terms  of  the  above  formulation, 
Pmax  of  the  TEMs  connected  in  parallel  has  the  same  characteristics 
as  the  series  connection,  as  long  as  the  conditions  of  Eqs.  (13)  and 
(16)  are  met.  In  Fig.  2,  the  AT  of  all  connected  TEMs  are  sorted  as 


Fig.  3.  Temperature  profile  of  the  hot  and  cold  fluids  with  parallel  heat  exchange  and 
the  TEG  critical  power  point. 


the  horizontal  axis,  whilst  Pmax  becomes  a  function  of  the  number 
of  connected  TEMs  corresponding  to  the  horizontal  position.  The 
TEMs  in  the  right  side  of  the  critical  (peak)  power  point  can  be 
regarded  merged  as  Vj  in  Fig.  1  as  well  as  (13)  and  (16)  according  to 
the  equivalent  circuit  theorem.  Obviously,  Pmax  does  not  mono- 
tonically  increase  with  the  increase  of  the  connected  TEMs,  and  in 
the  right  side  of  the  critical  power  point,  the  more  TEMs  with  small 
AT  are  connected  in  the  TEG,  the  less  power  can  be  obtained 
because  the  increase  of  the  equivalent  Rj  overwhelms  the  limited 
increase  of  Vj. 

For  a  TEG  with  the  aforementioned  quasi-one  1 D  heat  exchange 
configuration,  the  critical  power  point  can  be  defined  as  either  the 
position  of  a  TEM,  as  indicated  by  Lj  in  Fig.  3,  or  the  AT  of  the  TEM  at 
Lj.  From  0  to  Lj,  Th  keeps  decreasing  and  Tc  keeps  increasing  that 
after  Lj  the  temperature  difference  will  be  too  small  for  TEMs  to 
positively  generate.  If  the  fluid  length  is  fixed,  and  from  0  to  Lf,  Th 
keeps  decreasing  and  Tc  keeps  increasing,  the  temperature  differ¬ 
ence  distribution  between  them  will  be  shown  as  the  nearly  linear 
profiles  in  Fig.  4.  The  critical  power  point  can  be  defined  as  the 
slope  of  one  of  these  profiles,  where  the  TEM  at  If  has  zero 
contribution  to  the  TEG  power  output.  More  nonuniform  temper¬ 
ature  distributions  of  AT  will  make  one  or  more  TEMs  around  If  to 
become  depressed.  In  this  case,  the  critical  power  point  can  also  be 
defined  as  the  actual  value  of  A T(j)  at  the  position  Lf,  below  which 
the  temperature  difference  will  be  too  small  for  TEMs  to  positively 
generate.  For  the  2D  TEG,  the  critical  power  point  can  be  defined 
alike  as  the  AT  of  the  TEM  at  the  2D  position  (m  j)  in  Fig.  5,  where 
darker  bars  represent  depressed  TEMs  with  AT  coinciding  or 
smaller  than  the  critical  power  point.  Thus,  the  TEG  critical  power 
point  with  arbitrary  temperature  distribution  can  be  generalized  as 
the  AT  of  a  specific  TEM  position,  as  shown  in  Fig.  2.  When  the 
mismatch  power  loss  is  great  enough  to  enable  the  occurrence  of  a 
critical  power  point,  the  disconnection  of  one  TEM  or  a  certain 


Fig.  2.  General  characteristics  of  the  TEG  critical  power  point.  Fig.  4.  Temperature  difference  profiles  between  the  hot  and  cold  fluids  with  various 

ATi  >  AT2  >  ...  >  AT,  >  ...  >  ATm.  slopes  and  the  TEG  critical  power  point. 
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number  of  TEMs  with  A T  lower  than  the  critical  point  will  be 
favorable  to  the  overall  TEG  Pmax. 

3.  Experimental  diagnose  for  a  practical  TEG 

In  order  to  characterize  TEMs  in  various  interconnection  and 
temperature  configurations,  the  testing  apparatus  shown  in 
Fig.  6(al)  and  Fig.  7  is  able  to  independently  control  the  tempera¬ 
ture  difference  across  each  of  the  eight  TEMs  by  eight  independent 
electric  heaters.  All  the  data  were  obtained  using  TEM  devices  by 
Fuxin  Co.  (TEG1 -49-4.5-2.0-250).  Each  TEM  is  sandwiched  between 
a  hot  aluminum  heat  extender  and  the  cold  circulation  water 
chiller,  all  of  which  are  defined  as  a  TEM  unit.  In  every  TEM  unit, 
high  power  heating  elements  are  embedded  in  the  extender 
together  with  thermocouple  sensors,  and  in  this  way  precise  tem¬ 
perature  measurements  are  fed  back  to  temperature  controllers  to 


Fig.  7.  Illustration  of  the  TEG  with  setup  disassembly. 


automatically  fit  heater  temperatures  (7h)  towards  the  desired 
values  for  the  TEMs.  Thus,  the  temperature  mismatch  pattern  on 
the  TEM  array  can  be  arbitrarily  simulated.  Any  parallel/serial 
connection  of  the  TEMs  can  also  be  implemented  at  the  same  time 
by  combining  the  cell  minus/plus  lead  wires  but  using  the  identical 
stack  size. 

Given  the  good  flexibility  of  the  rig  to  enable  the  TEG  to  attain 
the  required  voltage  and  current  for  the  load,  the  2x4  paral¬ 
lel-serial  hybrid  TEM  matrix,  as  delineated  in  Fig.  6(a2),  is  actually 
connected  to  the  load  circuit  in  the  experiments.  A  minimal  boost 
DC-DC  converter  prototype  is  drawn  in  Fig.  6(a2)  to  represent  the 
power  electric  stage  part  in  the  load  circuit,  whilst  the  full  boost 
power  electric  schematics  actually  used  can  be  found  in  Ref.  27]. 
The  testing  results  of  the  hybrid  matrix  configuration  are  more 


Fig.  6.  Experimental  and  thermal-electric  coupled  modeling  architecture,  (al)  TEG  test  rig.  (a2)  Hierarchical  SPICE  modeling  at  the  system  level,  (bl)  TEMi  and  thermocouple 
closeup.  (b2)  Lumped  SPICE  modeling  at  the  TEM  level,  (cl)  Control  volume  of  the  p-type  and  n-type  thermoelements.  (c2)  SPICE  modeling  of  the  nonlinear  constitutive  equations 
of  thermoelectricity. 
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representative  and  generic  than  the  pure  parallel  and  serial  con¬ 
nections  of  TEM  arrays  because  the  latter  characteristics  are 
included  in  the  former  as  its  subset. 

The  accurate  characterization  of  the  TEG  critical  power  point  is 
firstly  obtained  by  sweeping  Thi,3,5,7  from  220  °C  down  to  80  °C, 
whilst  7h2,4,6,8  are  kept  constant  at  220  °C.  That  is,  the  initially 
uniform  Th  distribution  on  the  2x4  TEG  was  changed  to  be 
nonuniform,  where  one  half  TEMs  of  the  parallel  rows  are  subjected 
to  lower  temperature  than  the  other  half.  Owing  to  the  symmetric 
structure,  the  four  paralleled  TEM  pairs  in  the  series  circuit  are  of 
nearly  identical  thermal  conditions,  and  hence  the  power  loss 
within  the  mismatched  parallel  rows  is  mainly  examined. 

Fig.  8  shows  the  TEG  power  change  with  Thi,3,5,7  under  various 
resistive  loads.  The  output  of  the  TEG  can  also  be  connected  to  the 
electronic  load  through  a  DC-DC  converter  with  the  traditional 
MPPT  function  [27],  and  the  corresponding  results  are  presented  by 
the  curve  labeled  as  MPPT  mode,  where  the  equivalent  resistance  of 
the  load  circuit  automatically  coincides  the  TEG  internal  resistance 
by  adaptively  modulating  the  duty  cycle  of  the  DC-DC  converter. 

The  horizontals  in  Fig.  8,  labeled  with  “half’  and  corresponding 
to  various  load  modes  in  the  figure  legend,  represent  the  half  value 
of  the  initial  TEG  output  power  which  was  achieved  under  the 
maximum  Th  for  all  TEMs  (220  °C).  Any  point  of  intersection  be¬ 
tween  a  curve  and  its  counterpart  horizontal  with  the  same  load, 
marked  by  the  circles  in  Fig.  8,  represents  the  critical  power  point 
under  this  load.  The  value  of  Thi.3,5,7  at  this  point  indicates  that  AT 
for  TEMi^sj  has  been  small  enough  that  the  disconnection  of  them 
would  be  advantageous  to  the  TEG  power  because  TEM2,4,6,8 
themselves  along  are  now  able  to  produce  the  same  amount  of 
power  according  to  the  power  value  tested  under  the  initially 
uniform  Th.  If  Thi.3,5,7  further  decreases  that  AT  for  TEMi)3>5i7  is 
smaller  than  the  critical  power  point,  more  than  50%  of  original 
power  will  be  lost  even  half  TEMs  (TEM2,4,6,8)  are  still  provided  with 
the  high  temperature  220  °C. 

Without  changing  the  electrical  connection  inside  the  TEG,  the 
second  characterization  of  the  TEG  critical  power  point  is  obtained 
by  sweeping  Thi,2,3,4  from  220  °C  down  to  80  °C,  whilst  Th5,6,7,8  are 
kept  constant  at  220  °C.  In  this  case,  whether  it  is  a  high  temper¬ 
ature  or  a  low  temperature,  the  two  paralleled  TEMs  are  provided 
with  an  identical  Th  within  any  of  the  four  TEM  pairs.  Thus,  the 
mismatch  power  loss  of  the  TEM  string  is  mainly  examined. 
Compared  with  Fig.  8,  the  critical  power  points  under  different  load 
modes  in  Fig.  9  correspond  slightly  different  values  of  Thi.2,3,4,  but 
they  are  still  closely  around  90  °C.  In  both  Figs.  8  and  9,  the  output 
power  under  the  MPPT  mode  is  not  always  greater  than  other  fixed 
load  modes  due  to  the  limited  optimal  operation  range  of  the 
converter  [27],  i.e.,  the  power  consumption  of  the  converter  makes 


Fig.  8.  Output  power  and  critical  points  of  the  matrix  with  eight  TEMs  under  mis¬ 
matched  temperature  difference.  Load  condition  is  specified  in  legends. 


T hi,  T hi,  T h 3/  T hA  (degree) 

Fig.  9.  Output  power  and  critical  points  of  the  matrix  with  eight  TEMs  under  mis¬ 
matched  temperature  difference.  Load  condition  is  specified  in  legends. 


some  fixed  loads,  when  close  to  the  maximum  power  point, 
receiving  more  output  of  the  TEG  power  than  the  MPPT  mode.  The 
critical  power  point  characteristics  are  not  affected  under  the  MPPT 
mode  though,  and  it  is  validated  again  that  50%  power  is  lost  when 
half  TEMs  (TEMi)2,3,4)  are  provided  with  a  AT  equal  to  the  critical 
power  point. 

Although  the  eight  TEMs  composing  the  TEG  are  of  the  same 
model,  we  note  that  they  display  heterogeneous  electrical  perfor¬ 
mance  which  also  affects  the  mismatch  power  loss  and  the  critical 
power  point.  One  reason  is  that  these  TEMs  underwent  different 
generation  usages  (temperature,  current,  moisture,  and  compres¬ 
sive  load)  since  the  TEG  setup  had  been  built.  More  importantly,  the 
manufacturing  tolerance  of  low  cost  TEMs  such  as  TEG1 -49-4.5- 
2.0-250  can  lead  to  different  aging  rates  amongst  TEMs  [22].  The 
both  factors  are  typical  uncertainties  that  exist  in  practical  indus¬ 
trial  TEGs  nowadays,  which  make  the  accurate  prediction  by 
physical  simulation  difficult.  Including  the  uneven  aging  of  TEMs, 
however,  the  experimental  results  reported  above  still  clearly  verify 
the  theories  derived  in  Sec.2  that  the  critical  power  point  will  occur 
when  the  temperature  mismatch  arrives  at  a  certain  extent. 

4.  Generic  numerical  diagnose  by  SPICE  modeling 

With  the  rapid  development  of  thermoelectric  device  engi¬ 
neering  and  TEG  installation  standardization,  more  and  more  high 
end  TEMs  are  being  used  with  nearly  identical  aging  rates.  For  the 
new  generation  TEG,  an  accurate  and  systematic  TEG  model  is 
necessary  to  be  available  for  the  designer  to  determine  the  critical 
power  point  under  the  dynamically  changing  heat  source  temper¬ 
ature  distribution.  The  modeling  technique  based  upon  the 
fundamental  thermoelectric  physical  processes  can  also  provide  a 
reference  for  low  cost  TEGs  using  TEMs  with  heterogeneous  elec¬ 
trical  parameters  if  the  mismatch  power  loss  is  great  enough, 
despite  that  the  actual  critical  power  point  of  the  TEG  system 
consisting  of  TEMs  with  mismatched  properties  may  be  moderately 
different  from  the  modeling  result. 

The  theoretical  analyses  on  the  critical  power  point  in  Sec.2  are 
formulated  and  generalized  for  any  number  of  connected  TEMs.  It  is 
yet  mathematically  unrealistic  to  directly  apply  Eqs.  (13)  and  (16) 
due  to  the  nonlinearities  presented  within  Eqs.  (1),  (2),  (7)  and 
(8)  and  (5),  (6),  (9),  (10),  the  arbitrary  series-parallel-hybrid 
interconnection  of  practical  TEGs,  and  the  great  many  TEM  inter¬ 
connection  combinations  corresponding  to  various  critical  power 
point  definitions.  In  practical  control  of  a  partially  underheated 
TEM  array  caused  by  the  dynamically  changing  pattern  of  the 
temperature  profile,  both  the  fixed  TEMs  and  the  flexible  TEMs  for 
appropriate  removal/revival  are  multiple.  In  fact,  the  accurate 
performance  analysis  of  multiple  TEMs  has  been  identified  as  a 
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main  technical  issue  [21],  and  a  computer  model  is  indeed 
mandatory  to  numerically  determine  the  TEG  critical  power  point 
for  each  TEM  connection  combination. 

Moreover,  a  model-based  algorithm  is  also  important  in  adap¬ 
tively  perceiving  the  optimal  array  configuration  with  the 
maximum  output  power  in  terms  of  the  critical  power  point,  that  is, 
when  the  optimal  configuration  is  executed,  underheated  TEMs 
without  which  the  TEGS  can  output  more  power  can  be  automat¬ 
ically  removed,  and  removed  TEMs  can  also  be  automatically 
revived  if  their  previous  low  temperature  becomes  high  enough. 
Such  an  implementation  is  described  in  Ref.  [28]:  when  any  switch 
is  open,  its  relevant  TEM  will  be  removed  and  independent  from 
the  array,  whilst  the  removed  TEM  can  also  be  reconnected  into  the 
array  if  the  open  switch  becomes  closed.  When  the  temperature 
distribution  of  the  heat  source  is  nearly  uniform,  for  instance,  the 
modeless  trial  and  error  method  is  particularly  harmful  to  the 
generation  performance  as  the  continuous  trial  removal  of  certain 
TEMs  will  lower  the  output  even  the  removed  TEMs  can  be  revived 
later  on.  However,  the  undesirable  oscillation  and  loss  on  Pmax  can 
be  eliminated  by  the  model-based  algorithm  with  the  cost  of  using 
sensors. 

Therefore,  it  is  important  to  carefully  select  the  type  of  sensors 
in  terms  of  both  technical  and  economical  factors  in  order  to  fully 
take  the  advantages  of  the  model-based  algorithm.  Voltage  and 
current  sensors  can  be  used  in  terms  of  the  principles  formulated  in 
Sec.2  and  a  pure  circuit  network  model,  but  the  direct  measure¬ 
ment  of  V  and  R  of  a  TEM  in  the  array  is  inconvenient  and  may  cause 
oscillation  on  Pmax.  The  influences  of  the  thermal-electric  coupled 
multiphysics  effects  cannot  be  included  either  because  the  TEM 
currents  during  the  measurement  are  likely  different  from  the  real 
values  in  the  actually  operating  array,  leading  to  some  errors  of  V 
and  R  to  be  input  into  the  electrical  model  [29].  In  effect,  Eqs.  (1), 
(2),  (7)  and  (8)  and  (5),  (6),  (9),  (10)  clearly  indicate  that  V  and  R 
can  also  be  indirectly  computed  by  means  of  temperature  sensors 
only  through  an  electrothermal  modeling  approach  such  as  the  one 
in  Refs.  [24,25].  By  establishing  the  accurate  mapping  relation  be¬ 
tween  the  heat  source  temperature  distribution  and  the  TEG  power 
performance,  the  model-based  algorithm  can  still  decrease  the  cost 
brought  by  sophisticated  sensors  as  well  as  the  DAQ  (data  acqui¬ 
sition)  to  a  large  extent.  Table  1. 

Using  the  temperature  input  as  the  boundary  condition,  an  ideal 
TEG  model  to  be  used  in  the  determination  of  the  critical  power 
point  and  the  model-based  control  algorithm  should  have  the 
following  properties: 

1 )  The  model  is  based  on  the  physical  approach  and  the  universal 

thermoelectric  conversion  principles  rather  than  the  empirical 

approach  from  the  limited  existing  data. 


Table  1 

Definition  of  components  in  Fig.  6. 


Number 

Definition 

1-8 

TEM  units  1-8 

9 

Coolant  inlet 

10 

Coolant  outlet 

1.1 

Temperature  sensors  of  the  TEM  unit 

1.2 

Cartridge  heating  rods 

1.3 

Aluminum  heat  extender 

1.4 

TEM! 

1.4.1 

Top  ceramic  layer 

1.4.2 

Bottom  ceramic  layer 

1.4.3 

Thermoelements 

1.4.3. 1 

Conducting  bridges 

1. 4.3.2 

p-type  leg 

1. 4.3.3 

n-type  leg 
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2)  The  temperature  dependencies  of  ap<n,  pp>n,  and  /cp>n  of  the 
thermoelectric  materials  are  included. 

3)  The  multiphysics  coupled  effects  between  the  thermal  field  and 
the  electric  field  are  included. 

4)  The  ability  to  analyze  Pmax  of  various  series/parallel  in¬ 
terconnections  of  TEMs  with  different  temperature  boundaries 
(and  hence  different  V  and  R),  together  with  their  power  elec¬ 
tronics  stages. 

5)  A  user-friendly  design  environment. 

The  electrothermal  model  in  Ref.  [24]  is  based  on  the  full 
physical  details  and  nonlinear  constitutive  equations  of  the  ther¬ 
moelectric  conversion.  Fig.  6(bl)  illustrates  the  schematic  of  one 
TEM,  i.e.,  TEM1,  whose  fundamental  generation  thermocouple  is 
modeled  in  Fig.  6(b2),  where  the  potential  (VPin)  and  temperature 
(TPiI1)  distributions  inside  the  TEM  are  solved,  respectively.  Appro¬ 
priately  scaled  by  N,  Rc  [30]  is  simulated  in  the  electric  circuit,  and 
its  Joule  effect  is  divided  uniformly  to  be  simulated  in  the  hot  and 
cold  sides  of  the  thermal  circuit.  The  control  volume  of  the  ther¬ 
mocouple  is  depicted  in  Fig.  6(cl)  with  the  heat  and  electric  fluxes 
flowing  through,  whilst  its  SPICE  model  is  shown  in  Fig.  6(c2), 
where  both  the  temperature  dependent  heat  and  electricity  con¬ 
ductions,  as  well  as  the  multiphysics  coupled  effects,  are  numeri¬ 
cally  implemented.  Based  upon  electrothermal  analogy,  the  circuit 
simulator  SPICE  is  used  to  realize  the  multiphysics  coupled  algo¬ 
rithm,  where  the  thermal  and  electrical  simulations  are  performed 
in  the  same  modeling  tool.  As  shown  in  Fig.  6(b2)  and  (c2),  between 
the  isolated  thermal  and  electrical  circuits,  no  additional  routine  is 
required  to  control  their  information  exchange  and  iteration  for  the 
current-related  heat  source  terms,  i.e.,  Peltier,  Joule,  and  Thomson 
effects,  i.e., 

Sp,n(/)  =  +  V«p,n(rp,n(x))dxTp,n(x)/,  (17) 

*->p,n 

and  for  the  temperature  dependent  V  and  R.  By  the  TEM  model,  the 
corresponding  relationship  between  the  thermal/load  boundary 
condition  and  the  TEM  performance  parameters  (power,  efficiency, 
voltage  and  resistance)  has  been  elegantly  built  and  experimentally 
validated  [24]. 

By  upgrading  the  TEM  level  model  to  the  system  level  [25], 
multiple  TEMs  spatially  distributed  in  the  thermal  system  can 
freely  constitute  an  array  of  whatever  configuration,  including  both 
series  connection  and  parallel  connection.  The  inherent  nature  of 
the  SPICE  circuit  simulator  also  helps  the  implementation  of  the 
cosimulation  of  the  TEM  array  with  its  power  electronics  interface 
and  MPPT  controller.  In  the  TEG  model,  a  hierarchical  and  TEM- 
object  oriented  strategy  is  utilized  that  each  TEM  type  is  defined 
as  a  class,  in  which  all  detailed  simulation  procedures  and  device 
data  of  the  TEM  type  are  packaged  and  become  invisible  to  the 
outside  of  the  class.  At  the  system  level,  only  one  TEM  class  is 
invoked  for  all  TEM  instances  of  the  same  type,  and  every  TEM 
instance  only  has  four  input  and  output  (I/O)  ports  to  be  manipu¬ 
lated,  i.e.,  Th,  Tc,  plus,  and  minus  lead  wires,  as  shown  in  Fig.  6(a2). 
In  this  way,  the  maintainability,  scalability,  development  efficiency, 
and  execution  efficiency  of  the  TEG  modeling  code  is  essentially 
improved  [25]. 

Using  the  SPICE  modeling,  the  TEG  setup  described  in  Sec.3  and 
Fig.  6(a)  is  simulated  to  numerically  quantify  the  characteristics  of 
the  critical  power  point.  Firstly  a  DC  sweep  analysis  on  the  voltage 
source  emulating  T m  is  performed  to  obtain  the  critical  power  point 
of  a  single  low  temperature  TEM.  In  Fig.  10,  it  is  noteworthy  that  the 
output  power  is  not  completely  constant  when  TEM4  is  excluded 
from  the  SPICE  simulation  because  the  water  temperature  variation 
is  taken  into  account  at  the  cold  side.  As  shown  in  Fig.  6(a2),  with 
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0.7  Ohm 


T  h4  (degree) 

Fig.  10.  Output  power  and  critical  points  of  the  matrix  with  TEM4  under  low  Th. 
Thi— 3,5— 8  =  220  °C. 


0.7  Ohm 


T  h4(degree) 

Fig.  12.  Output  power  and  critical  points  of  the  matrix  with  TEM4j8  under  low  Th- 
Thi-3,5-7  =  220  °C.  Th8  =  60  °C. 


the  decrease  of  Th4,  the  heat  flow  added  from  TEM4  to  the  coolant 
fluid  is  decreased,  too.  Consequently,  Tc  for  the  downstream  TEMs  is 
gradually  lowered  to  furnish  these  TEMs  larger  AT,  and  the  TEG 
output  power  is  hence  increased. 

Fig.  10  shows  that  the  critical  power  point  value  of  Th4  is  about 
80  °C  and  60  °C  for  R\  =  1.5Q  and  R\  =  1.1  Q,  respectively.  When 
R\  =  0.7Q,  the  critical  power  point  does  not  appear  within  the  figure 
scope,  which  indicates  that  AT  for  TEM4  is  still  not  small  enough 
and  the  disconnection  of  it  would  decrease  P.  In  Fig.  11,  Th  of  other 
TEMs  is  purposely  designed  to  be  nonuniform  to  enhance  the 
testing  diversity.  However,  the  critical  power  point  characteristics 
are  nearly  identical  to  Fig.  10  in  this  case. 

The  second  DC  sweep  analysis  is  performed  to  obtain  the  critical 
power  point  characteristics  of  two  low  temperature  TEMs,  i.e., 
TEM4  and  TEM8.  In  Fig.  12,  Ths  is  fixed  at  60  °C,  and  Th4  is  decreased 
from  220  °C  down  to  60  °C.  It  is  found  that  the  critical  power  point 
of  Th4  is  about  120  °C  and  80  °C  for  R\  =  1.5Q  and  R\  =  1.1  Q, 
respectively,  and  when  Pi  =  0.7Q,  the  critical  power  point  does  not 
appear  within  the  figure  scope.  When  both  Th4  and  Ths  are  swept 
from  the  high  temperature  to  the  low  temperature,  the  critical 
power  point  is  about  75  °C  and  60  °C  for  Pi  =  1.5Q  and  R\  =  1.1  Q,  and 
again,  the  critical  power  point  does  not  appear  when  R\  =  0.7Q,  as 
shown  in  Fig.  13. 

5.  Conclusions 

The  optimization  approach  in  terms  of  the  critical  power  points 
is  an  extremely  important  research  topic  in  the  development  of  TEG 
applications.  Such  an  optimization  can  reduce  the  negative  effects 
of  the  temperature  mismatch  on  the  output  power,  which  has  been 
commonly  observed  by  TEG  system  designers,  and  hence  consid¬ 
erably  increase  the  actual  energy  efficiency  of  practical  TEGs.  One 


key  step  in  making  the  optimization  feasible  is  to  determine  the 
critical  power  points  and  quantify  their  magnitude  under  various 
operating  conditions.  Exploiting  the  principles  of  battery 
mismatch,  rigorous  analytical  derivation  is  given  in  this  paper  on 
the  critical  power  point.  It  is  mathematically  proven  that  when 
some  TEMs  are  exposed  to  a  smaller  heat  source  AT  in  a  TEM  array, 
the  generated  TEG  power  can  be  more  by  removing  them  because 
of  their  internal  energy  consumption.  The  theoretical  conceptions 
and  conditions  of  the  critical  power  point  are  illustrated  by  ex¬ 
periments  and  numerical  simulations  of  a  series-parallel  hybrid 
interconnected  TEM  matrix,  respectively.  In  experimental  and  nu¬ 
merical  results,  a  number  of  critical  power  points  are  disclosed  for 
the  TEG  prototype,  where  Th  mostly  ranges  from  120  °C  to  60  °C. 
The  experimental  results  prove  themselves  that  the  critical  power 
points  for  the  TEG  composed  of  TEG1 -49-4.5-2.0-250  are  accurately 
determined,  and  critical  power  points  for  other  systems  can  be 
accurately  determined  by  the  same  approach.  More  generally,  the 
SPICE  simulator  is  used  as  the  modeling  and  programming  tool  and 
facilitates  the  maximum-power  oriented  configuration  for  live 
TEGs  at  large  by  means  of  the  hierarchical  cosimulation  approach, 
where  all  major  thermoelectric  nonlinearities,  i.e.,  thermo-electro 
coupled  multiphysics  behaviors  and  temperature  dependent  ma¬ 
terial  properties,  are  included.  The  accuracy  of  the  TEG  model  has 
also  been  verified  by  many  previous  experiments  [24-26,29], 
including  those  on  the  critical  power  points  [28]. 

In  principle,  it  is  possible  that  removing  certain  TEM(s)  of  higher 
temperature  difference  can  give  more  output  power  than  removing 
certain  TEM(s)  of  lower  temperature  difference  from  the  same 
configuration.  Such  examples  are  particularly  likely  to  be  illustrated 
when  the  TEM  array  is  of  asymmetric  and  irregular  structure  and/or 
when  the  TEM  array  is  operating  in  connection  with  other  power 


T  h4  (degree) 
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Fig.  11.  Output  power  and  critical  points  of  the  matrix  with  TEM4  under  low  Th.  Fig.  13.  Output  power  and  critical  points  of  the  matrix  with  TEM4  8  under  low  Th. 
Tm-3,5-8  =  220  °C,210  °C,200  °C,190  °C,180  °C,170  °C,160  °C.  Thi_3i5_7  -  220  °C,210  °C,200  °C,180  °C,170  °C,160  °C. 
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sources  such  as  fuel  cells/PV  cells/battery  packs,  where,  apart  from 
the  temperature,  the  position  of  each  TEM  also  matters.  In  this  case, 
an  exhaustive  traveling  of  all  possibilities  of  power  increase  is  the 
safest  way  to  attain  the  globally  optimal  configuration.  It  is  yet 
important  to  acknowledge  that  the  conception  of  critical  power 
point,  to  search  only  within  those  topologies  removing  under¬ 
heated  TEMs  instead  of  all  topologies,  is  valuable  in  dealing  with 
common  TEGs  which  are  mostly  applied  in  practice.  The  possibility 
will  also  exist  to  work  on  a  range  of  problems  relating  the  critical 
power  point  to  TEG  security,  as  the  physically  destroyed  TEMs  may 
also  be  adaptively  removed  so  that  the  performance  of  the  TEG  is 
actually  able  to  self-heal. 
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Nomenclature 

S  heat  source  term  (W) 

Dw  area  of  the  heat  source  (m2) 

dx  ID  control  volume  length  in  thermoelements  (m) 

I  current  (A) 

i  TEM  number 

k  thermal  conductivities  (W/mK) 

L  length  of  the  p-  and  n-type  legs  (m) 

If  length  of  the  fluid  path  (m) 

Lj  position  of  the  TEM  at  the  critical  power  point  (m) 

M  total  number  of  TEMs  in  the  TEG 

m  maximum  number  of  TEMs  in  the  horizontal  coordinate 

of  2D  TEG 

N  number  of  pairs  of  thermoelements  in  each  TEM 

n  maximum  number  of  TEMs  in  the  vertical  coordinate  of 

2D  TEG 

P  power  on  the  load  (W) 

Pmax  maximum  power  (W) 

Q.  conduction  heat  flow  in  p-type  or  n-type  thermoelement 

(W) 

R  internal  resistance  (Q) 

Rc  contact  resistance  across  conducting  bridges 

R\  load  resistance  (Q) 

S  uniform  cross-sectional  area  of  TEM  (m2) 

Sn  uniform  cross-sectional  area  of  n-type  leg  (m2) 

Sp  uniform  cross-sectional  area  of  p-type  leg  (m2) 

T  temperature  (K) 

V  open-circuit  electric  potential  (V) 

w  dimension  of  each  TEM  in  the  direction  of  the  fluid  flow 

(m) 

Greek  symbols 

a  Seebeck  coefficient,  (VK-1) 

AT  temperature  difference  across  the  TEM,  Th  -  Tc,  (K) 
p  electrical  resistivity  (Q  m) 

Subscripts 
c  cold  side 

h  hot  side 

j  TEM  number 

n  n-type  material  or  leg 

p  p-type  material  or  leg 
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